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^ ■ Abstract 
O 

OS , We apply a generalized numerical transfer-matrix method to the two- 



dimensional Ising ferromagnet in a nonzero field to obtain complex con- 
strained free energies. Below T c certain eigenstates of the transfer matrix 



C. 

O , are identified as representing a metastable phase. The imaginary parts of the 

O 

metastable constrained free energies are found to agree with a field-theoretical 
droplet model for a wide range of fields, allowing us to numerically estimate 
the average free-energy cost of a critical cluster. We find excellent agreement 
with the equilibrium cluster free energy obtained by a Wulff construction with 
the exact, anisotropic zero-field surface tension, and we present strong evi- 
dence for Goldstone modes on the critical cluster surface. Our results are 
also fully consistent with average metastable lifetimes from previous Monte 
Carlo simulations. The study indicates that our constrained-transfer-matrix 
technique provides a nonperturbative numerical method to obtain an ana- 
lytic continuation of the free energy around the essential singularity at the 
first-order transition. 
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Despite the fact that metastability is commonly observed in systems ranging from flu- 
ids, magnets, and alloys, to the "false vacuum" in high-energy physics, a fully satisfactory 
microscopic description remains elusive. Although metastable states do not globally mini- 
mize the free energy, they persist for periods many orders of magnitude longer than other 
characteristic timescales of the system jrj,|2[ . This observation has prompted efforts to study 
the statistical mechanics of metastable states as if they were equilibrium phases, using con- 
strained partition functions that exclude or severely reduce the probabilities of microstates 
that dominate in equilibrium Q. By applying such ideas to a field-theoretical droplet model 
with Fokker-Planck dynamics, it has been shown [|3],f|] that the nucleation rate for droplets of 
the equilibrium phase in a metastable background is proportional to the imaginary part of a 
complex-valued constrained free energy obtained by analytic continuation from the equilib- 
rium phase into the metastable phase. This result was derived for ultra- weak fields only [|J, 
and despite extensive subsequent work [|5|-[7| its domain of validity is not completely clear. 

Recently one of us introduced a numerical technique, the Constrained- Transfer-Matrix 
(CTM) method, which extends transfer-matrix calculations to metastable and unstable 
states H , providing a means to compute complex- valued "constrained free energies" . How- 
ever, in contrast to the field-theoretical droplet model the CTM method does not 
necessitate the explicit introduction of droplets. In this Letter we apply this method to the 
two-dimensional square-lattice nearest-neighbor Ising ferromagnet in a magnetic field and 
below its critical temperature. Metastable states in short-range force systems have finite, al- 
beit very long, lifetimes |6|j9|-[T2| , and it is now quite well established that the free energy has 



an essential singularity at H=0 fl^]T3}j . This contrasts sharply with the infinitely long-lived 
metastability displayed by systems with weak long-range interactions in the thermodynamic 
limit p|,p|,|7|,p| fl5| , p!6[] and the corresponding spinodal branch-point singularity. The numeri- 
cal results for the metastable free energy that we present here are in excellent agreement 



with series expansions [|T7|-|20|j and Monte Carlo simulations |IT|JT^ , |2T|J2^| , and indicate that 
the CTM technique provides a nonperturbative numerical method to obtain an analytic 
continuation of the free energy around the essential singularity at H=0. 



The Ising model is defined by the Hamiltonian 7i = —JJ2(i,j)° r i°'j ~ HJ2i°~i, where 
(Tj = ±1 is the spin at site i, J>0 is the nearest-neighbor coupling, the first sum is over all 
nearest-neighbor pairs on the square lattice, and H is the field. 

The conventional transfer-matrix formalism |23| provides a method to calculate the equi- 
librium partition function and correlation lengths for an Nxoo system from the eigenvalues 
of the transfer matrix To- This matrix is defined by dividing the system into one-dimensional 
layers of N spins each and by decomposing the total Hamiltonian (in this case TC from above) 
into partial Hamiltonians H(Xi,Yi + i), giving (J^"/|T |V/+i) = exp[— ^H(Xi, Vj+i)]. Here T 
is the temperature in energy units, and Xi and Y/ + i are the configurations of the l-th and 
(/ + l)-th layer, respectively. Since To is a positive matrix, its dominant eigenvalue is posi- 
tive and nondegenerate (i.e., A > \X a \ for all a>0), and the corresponding eigenvector |0) 
is the only one whose elements are all positive p3| . The CTM method extends the scope 



of transfer-matrix calculations to metastable and unstable states by defining constrained 
probability densities constructed from the eigenvectors \a) of T , 

P a {X h Y m ) = {a\X{){X l \{\~ 1 Tj m - l \\Y m ){Y m \a) 

P a (X t ) = (alX^X^a) . (1) 

The matrix Tq, is chosen to commute with To, but has X a as its dominant eigenvalue, 
rather than A . This ensures that limi^i^oo P a (X h Y m ) = P a (Xi)P a (Y m ), corresponding to 
stochastic independence of widely separated layers. As a consequence, the T a for a>0 are 
not positive matrices. As in Refs. |8|,^,p^l we use 

T a = A J £ \(3)^((3\ + \a)(a\ + £ |/?>^</?i) . (2) 

From T Q one obtains constrained free-energy densities, 

F « = 4 E P„(Xt,Y l+1 )H(X u Y l+1 ) -TS a , (3a) 

Xi,Yl+i 



where 



S « = -7j E ^(^,^+i)Ln(^|A^ 1 T a |^ +1 ) (3b) 

Xi,Y l+1 



are constrained entropy densities, defined in analogy with the source entropy of a stationary 
Markov information source |^3J . Negative elements of T a give rise to the imaginary part of 



F a through Lnz, the principal branch of the complex logarithm. The quantities defined in 
Eqs. (|2])-(]3"b"|) reduce to their standard equilibrium values |23|] for a—0. 



To calculate F a from Eq. ( |3"a| ) one has to sum over all the eigenstates of T . To accomplish 
this with sufficient accuracy for large N, diagonalization routines in 128-bit precision for the 
Cray Y-MP/432 vector-supercomputer, involving block-diagonalization of To, were custom- 
written. For T/J< 0.4, the attainable values of N were limited by numerical underflow 
for small \H\, whereas for higher T they were limited by the available memory. The total 
computer time spent on this study was on the order of 800 CPU hours. 

Following previous transfer-matrix studies f5|||], the metastable state for any particular 
H>0 was identified as the eigenstate \a) whose magnetization is closest to —1.0. Fig. 1 shows 
lmF a vs. J/H for those a that contribute to the metastable state for N=9 and 10. Each of 
the lobes corresponds to a different a. Their minima, especially for small H, have extremely 
small values. The crossings correspond to avoided crossings in the eigenvalue spectrum of T . 
At T=0 the exact values of H at which these crossings occur are H=0 and H/ J = 2/(N—m) 
with m— 1, . . . , N—l. For the temperatures studied, the field values of the crossings deviate 
by at most 3% from their values at T=0. As pointed out previously |5|,§], these crossings 
should correspond to zeros of the constrained metastable partition function, and as N^oo 
the diverging density of crossings near H=0 gives rise to an essential singularity [ fPfj . 



We find two different scaling regimes for the minima of the metastable lmF a . The 
minimum which lies in the field interval < H/J < 2/(N— 1) decays exponentially with 
N. The value of H at which that minimum occurs, is proportional to N~ 2 , whereas the 
corresponding dominant length scale is proportional to N. For H/ J between approximately 
2/(N — 2) and 2, the metastable ImF a , while still extremely small, are independent of N. 
These results can be explained by a droplet model in which the initial step in the decay of 
the metastable state is the random formation of clusters that exceed a critical size. 

In weak fields, the average critical cluster is larger than N. Therefore, the free-energy 



cost of creating a critical cluster is simply the free energy of an interface cutting across 
the system, and thus proportional to N. Assuming that the Boltzmann factor of this free- 
energy cost, as well as the minimum of the metastable lmF a , is inversely proportional to the 
nucleation rate, one finds that lmF a decays exponentially with N. The N~ 2 dependence of 
the field at which this minimum occurs can also be derived by a droplet-model argument 



Field-theoretical droplet-model calculations for infinitely large systems in the limit of 
very weak magnetic fields |p|,p6| show that in two dimensions the imaginary part of the 



metastable free energy should be given as [20 



ImFrr 



1 



B\H\°exp 



T 4 Iff AM 



0(ff 2 



(4) 



where B is a nonuniversal constant. The leading term in the exponential is 1/T times the 
nonuniversal free-energy cost of a critical cluster, with £ 2 being the square of its surface free 
energy divided by its area, and AM the difference in magnetization between the equilibrium 
and metastable states. The prefactor exponent 6 is determined by Goldstone excitations on 
the critical cluster surface and is expected to be universal. For d=2, 6 is expected to be 
unity if Goldstone modes are present JTT| . p~5| . pD|J^B| , p7|] and —1 if they are absent [|T8(] . Series 



expansions fli7H19|1 and Monte Carlo simulations pl| support Eq. dJ) with 6=1. Furthermore, 
there is substantial evidence |]20|,|2l| that, in the limit of weak fields, S 2 is given by its 



equilibrium value, which may be obtained by combining a Wulff construction with the exact, 
anisotropic zero- field surface tension |28j , and that AM is given by twice the exact zero-field 



magnetization 20 



Despite the fact that Eq. (£|) was derived in the weak- field limit, Fig. 1 indicates that 
it also gives the correct leading-order behavior of the minima of the metastable lmF a for 
intermediate fields, 2/ (N— 2) < H/J<2. The dashed straight line in Fig. 1 is drawn through 
the two minima between J/H=3 and 4 for iV=10 only. However, it also connects several 
minima at smaller values of J/H for both iV=10 and 9, and its slope gives a rough estimate 
for S 2 /AM in Eq. (|j). In this field region the size of the average critical cluster is smaller 



than the finite width of the system. Indeed, we find that the average droplet diameters, 
as determined from the droplet model, range from about 8 lattice units for N=10 near the 



cross-over field if/ J ~ 2/(10— 2), down to 1 near H/ Jm2 for all T studied p5| . Furthermore, 
for all temperatures and all fields H/ J<1.2, we find that both the metastable and the stable 
single-phase correlation lengths, as determined from the transfer-matrix level spacings, are 
smaller than the average droplet diameters. Except for the largest fields studied, this ensures 
sharp cluster surfaces. 

In order to determine E 2 /AM and b in a more systematic way, we performed a three- 
parameter linear least-squares fit of the logarithm of the minima of the metastable lmF a 
to the form of Eq. ([|). It must be emphasized, however, that since our results for lmF a 
are numerically exact, the least-squares method was used merely as a tool to estimate the 
parameters numerically. We have therefore assigned error bars such that \ 2 P er degree of 
freedom is approximately unity. 

Our results for E 2 /AM are in agreement with the hypothesis that £ 2 is equal to its 
equilibrium value obtained by combining a Wulff construction with the exact, anisotropic 
zero-field surface tension P5[ , and that AM equals twice the exact zero-field equilibrium 
magnetization. We find agreement to within 10% for all temperatures. We therefore fix 
E 2 /AM to its equilibrium value and obtain b from a two-parameter linear least squares fit. 
For T j J>0.4, the results are consistent with the presence and inconsistent with the absence 
of Goldstone modes on the surface of the critical clusters, with a value of 6=0.87 ± 0.19 at 
T j J= 1.0. The limits on the attainable values of N prevent us from probing sufficiently small 
fields to obtain reliable estimates of b for T/J=0A. Also, large finite-size effects prevent us 
from attaining sufficiently high temperatures to verify any effects of the deviation of AM 
from 2 [0. 

Fig. 2 shows a comparison between the equilibrium value of £ 2 /AM with its value 
obtained from a two-parameter fit, in which b is fixed to 1. It strongly supports the idea 
that the metastable state decays through the formation of critical clusters. These are on the 
average almost square at low temperatures and nearly circular at higher temperatures. 



Our results are also consistent with direct measurements of metastable lifetimes by Monte 
Carlo simulation P,p|-p^,p^] , in which three different scaling regimes are found. The two 
low-field regimes correspond directly to those found here for the metastable IulFq,. 

In summary, the Constrained- Transfer-Matrix method || applied to the two-dimensional 
square-lattice Ising ferromagnet allows us to numerically obtain complex metastable free 
energies. Their imaginary parts depend on system size and magnetic field in a way fully 
consistent with the interpretation that this technique provides a nonperturbative numerical 
method to obtain an analytic continuation of the free energy from the equilibrium phase 
into the metastable state around the essential singularity at H=0. However, the CTM 
method does not rely on the explicit introduction of clusters of the stable phase in the 
metastable background. Comparing IhlF q with field-theoretical droplet-model calculations 
of the same quantity H J20|J26[| , we estimated the average free-energy cost of a critical cluster 
over a wide range of fields and temperatures, finding excellent quantitative agreement with 
theoretical predictions |2j],|2lJ] involving the surface free energy of equilibrium cluster shapes 
28fl . The above comparison also gives strong evidence for Goldstone modes on the critical 



cluster surfaces [^,|I^,^,^T|,^,^] . A more detailed account of this study will be reported 
elsewhere 11251 • 
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FIGURES 

FIG. 1. The imaginary part of the metastable constrained free energy, ImF a , vs. J/H for 
two different system sizes at T/J=1.0 (T/T c ~0.441). The metastable Imi^, which is expected 
to be proportional to the nucleation rate, varies over twelve orders of magnitude. There are two 
different scaling regions. For small \H\, ImF a vanishes exponentially with N, and for large H it is 
independent of N. The straight dashed line is drawn through the two minima marked by stars in 
the figure for N=W only, but it also connects minima at smaller values of J/H, both for iV=10 
and iV=9, as expected from Eq. (Q). 

FIG. 2. Plot of X 2 /AM vs. T/J, as obtained from a two-parameter linear least-squares fit of 
the logarithm of ImF a to Eq. @, setting b=l. The solid line corresponds to the equilibrium surface 
free energy calculated using a Wulff construction with the exact, anisotropic, zero- field Onsager 
surface tension cr(T) [^] divided by twice the exact zero- field magnetization M(T). The cluster 
shape given by £ 2 interpolates between a square at T=0, given by 8a 2 /M, and a circle at T=T C , 
given by 27r<7 2 /M, both shown for the whole temperature range as dashed lines. The metastable 
transfer- matrix results follow the equilibrium curve closely. The dotted vertical line marks T c . 
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